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An ion beam analysis system was established on a 1.7 MV tandem accelerator, enabling Rutherford backscat- 
tering (RBS), elastic recoil detection (ERD), nuclear reaction analysis (NRA) and channeling measure- 
ments. The system was tested by performing qualitative and quantitative analysis of Si, Ni/Si, BiFeO,:La/Si, 
MoC/Mo/Si and TiBN/Si samples. RBS of a BiFeO,:La film was used as system calibration. Tested by ion 
beam channeling, a Si(100) is of good crystallinity (Xmin = 3.01%). For thin film samples, the measured thick- 
ness agrees well with simulation results by SIMNRA. In particular, composition of a MoC/Mo/Si and TiBN film 
samples were analyzed by RBS and non-Rutherford elastic backscattering. 
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I. INTRODUCTION 


Based on interactions of energetic ions with target materi- 
als, ion beam analysis techniques are widely used in charac- 
terization of thin films and multilayer structures. They pro- 
vide quantitative information of elemental components of the 
solid from several atomic layers to micrometer depth from the 
surface, or depth profiling of impurities. Therefore, it attract- 
ed much attention in the ion beam community and many ion 
beam systems have been developed and upgraded. These in- 
clude conventional ion beam analysis, microbeams [1—4], and 
combination of ion beam systems, such as the 5 MV tandem 
accelerator interfaced with the 200 MeV proton synchrotron 
in the Wakasa Wan Energy Research Center [5] and the 
multi-ion irradiation platform JANNUS (Joint Accelerators 
for Nano Sciences and Nuclear Simulation) at Saclay where 
the first triple beam irradiation using Fe, He and H ion beams 
was carried out in 2010 [6]. 

Ion beam analysis methods are frequently used due to 
their efficiency in providing a broad spectrum of qualitative 
and quantitative information on elemental composition and 
crystal structure of solids. While Rutherford backscatter- 
ing (RBS) provides composition and depth profiles of heavy 
elements in a light matrix, elastic recoil detection analysis 
(ERDA) is unique in profiling light elements (C, B, Li, H 
etc.) in a heavy matrix. Nuclear reaction analysis (NRA) 
provides specific elemental profiling which can resolve iso- 
topes. Proton induced X-ray emission (PIXE) enables simul- 
taneous detection of elements of intermediate atomic num- 
bers. Ion beam channeling, compared with random ion beam 
scattering spectroscopy, gives high precision measurement 
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of crystallinity of single-crystalline wafers or epitaxial lay- 
ers [7]. Ion beam analysis methods have a common advan- 
tage of non-destructiveness in the sense that the correspond- 
ing ion-solid interactions do not consume the target material, 
since the beam current is at analytical level, i.e., a few nA. 
Damage may result only, for sensitive samples, from ioniza- 
tion or heating by the ion beams, sputtering or nuclear stop- 
ping and in the case of channeling studies in single crystals, 
from defect production inducing dechanneling [8]. 

The accelerator system in Wuhan University is multifunc- 
tion equipment, namely, the 1.7 MV tandem accelerator is in- 
terfaced with a 200 kV ion implanter and a transmission elec- 
tron microscope (TEM). The ability of ion beam technique is 
extended to emerging fields of physics and materials science. 
In the present work, we have designed RBS/PIXE/C system 
aiming at characterization of solids with a broad range of el- 
ements in the same chamber and in situ measurement of the 
implanted sample in a dual beam chamber. The advantage 
of low-energy RBS technique lies in its capability of quan- 
titative analysis of major and minor constituents lying in the 
first 0.5-2.0 um of a material depending on the sample struc- 
ture and composition. And the detection limits for some light 
elements may further be improved by using resonant scatter- 
ing. In this paper we report the tests of basic functions of the 
system based on RBS and channeling, and demonstrate quan- 
titative characterization of elemental profiling, and interface 
analysis, and crystal quality of wafers and thin film materials. 


Il. SYSTEM DESIGN 


Figure 1 shows a schematic diagram of the in situ ion beam 
system. Ion beams extracted from the negative ion source by 
cesium sputtering (SNICS) are transported to the low energy 
cluster ion beam chamber for low energy ion implantation 
or to the RBS/PIXE/C chamber after transportation lines of 
the tandem accelerator. The dual beam chamber is used for 
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in situ implantation with ions from the tandetron or 200 kV 
implanter. It will also be possible to conduct in situ ion beam 
analysis of the implanted samples with the establishment of 
the RBS/PIXE/C system in the dual beam chamber. The ion 
beam analysis line thus consists of the SNICS, tandetron and 
analysis chamber. 
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Fig. 2. (Color online) Schematic diagram of the measurement sys- 
tem of the ion beam analysis line. 


A schematic diagram of the measurement system is shown 
in Fig. 2. Conventional ion beam scattering, recoil detection 
and nuclear reaction are performed with Si:Au surface barrier 
detectors (SBD) connected through a charge sensitive pream- 
plifier, an amplifier and a multichannel analyzer (MCA). The 
total number of ions is recorded by a charge integrator. For 
channeling measurement, a single channel analyzer (SCA) 
and a linear rate-meter are used. For PIXE measurement, we 
have an Ortec SLP-1080 detector with liquid nitrogen cool- 
ing and a miniature detector XR-100CR with thermal elec- 
tric cooling. The SLP-1080 detector is put into the chamber 
through a cinlyndrical port. The XR-100CR detector is placed 
in the chamber, and the electronics for the RBS measurement 
can be used for signal amplification and recording. 

Ion beams extracted from the SNICS are focused and de- 
flected before they were accelerated to 1-3 MeV, depending 
on energy requirement for different measurements. The ac- 
celerated ion beams are focused again by the electrostatic 
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quadrupole lens and deflected into the analysis chamber by 
a magnetic switcher. The samples are placed on a 670mm 
copper holder mounted on a four-axis goniometer, which 
is computer controlled through a stepping motor controller. 
During the measurement, the stepping motor controller con- 
trols movement of the goniometer, providing three orthogonal 
rotations (two tilts and one azimuth) and one translation (y- 
axis). The z-axis is aligned with the incident ion beam direc- 
tion. The y-axis is defined by orthogonality with both the x- 
and z-axis, which complete a right-hand triad. Specification 
of the sample orientation is described by three angles, y, n 
and C, corresponding to rotation of the x-, y- and z-axes, re- 
spectively. The goniometer enables the samples to move a- 
long the y-axis or rotate around the x-, y- and z-axes. The s- 
patial relationship between the sample and incident ion beam 
is described by the angle 6 = (7° + ¢7)'/? and the azimuthal 
angle @. 

In RBS analysis, the ion beam was collimated to 0.3mm 
size by an aperture slot. The SBD, with an effective area of 
50mm”, is positioned at a distance of 10-15 cm to the sam- 
ple and 150°—170° to the beam incidence. The signals, af- 
ter preamplifier and the spectroscopic amplifier, were sent to 
the MCA installed in the computer. The current integrator 
records the charge collected during the experiment. 

The RBS/C measurement has three more steps. First, sig- 
nals from the amplifier are sent to a gate interface of the MCA 
through a single channel analyzer (SCA), simultaneously, an 
energy window is selected through coincidence by regulating 
the upper and lower level discriminators of the SCA. Next, 
the gate interface is disconnected and signals of the select- 
ed energy range from the SCA are sent to the linear rate- 
meters. Finally, tilt or rotation of the crystal sample relative to 
the beam incidence is controlled by the stepping motor con- 
troller, and the reduction of scattering yield is displayed on 
the linear rate-meters to decide channeling effect of the inci- 
dent beams [9]. Detailed information of how to determine the 
channel of crystal sample by computerized control program 
was described in the literature [10]. 


Il. ANALYSIS AND DISCUSSION 
A. System calibration and determination of crystal quality 


The RBS system was calibrated by measuring the Au peak 
position and Si edge of a Au-deposited Si sample. RBS 
measurement of a BiFeO,:La layer sample was performed. 
BiFeO, is a typical multiferroic material with ferroelectric 
and antiferromagnetic properties, has and is of broad appli- 
cation prospects in spintronics and multibit memory devices. 
Doping it with rare earth elements may significantly improve 
its ferroelectric and ferromagnetic properties. It is important 
to measure each element accurately. A thin film of La-doped 
BiFeO, was prepared by sol gel method on a Si(100) sub- 
strate. Fig. 3(a) shows RBS spectrum of the sample using 
2.9MeV ‘Li** beams with a backscattering angle of 160°. 
The leading edge of Si is at channel 147 (at half-height) and 
the Bi, La, Fe and O peak positions are at channel numbers of 
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Fig. 3. (Color online) RBS spectrum of La-doped BiFeO, film on Si 
(100) (a) and RBS/C spectrum of single-crystal Si (100) (b). 


386, 362, 265, and 62, respectively. They correspond to the 
kinematic factors of 0.8778 for Bi, 0.8219 for La, 0.6126 for 
Fe, 0.3713 for Si, and 0.1607 for O. This indicates that Li ion 
beam has excellent mass resolution for many elements of film 
analysis. Assuming that the deposited La-BiFeO, thin layer 
is distributed uniformly on the substrate, the film thickness is 
estimated at 135 nm. Compared with other methods, the RBS 
analysis is very accurate, with an error of just 2%. 


Figure 3(b) shows RBS/C measurement results of a sin- 
gle crystal Si(100) sample using 2.75 MeV Li ion beams. 
Compare to the random spectrum, the backscattering yield 
greatly reduced in the aligned spectrum, the flat shape of 
which demonstrates that the incident ions go deeply into the 
crystal channels without dechanneling. At the high-energy 
edge, a small peak was observed, revealing an amorphous lay- 
er in the surface owning to the cutting during sample prepara- 
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tion. The minimum yield, ymin, defined as the backscattered 
yield ratio of the aligned spectrum to the random spectrum 
in the same channel region, i.e., the same depth, represents 
structural quality of the whole lattice of a monatomic crystal. 
For a perfect crystal, Xmin could be 1-2% [9]. The measured 
Xmin for the Si(100) sample is 3.01%, indicating a high quality 
crystal structure. 


B. Quantitative analysis of thin films 


Thickness of a Ni film prepared by electron beam evapora- 
tion on SiO,:Si substrate was measured by RBS with 2.9 MeV 
TLi?* in normal incidence at a backscattering angle of 160° 
(Fig. 4). The yield of backscattered Li ions was detected 
by a Si:Au surface barrier detector with an effective detec- 
tion area of 50 mm?, positioned at 112.6 mm from the sam- 
ple, with a solid angle of Q = 3.94 msr. The measured spec- 
trum overlaps with the theoretical spectrum simulated using 
the iterative analytical code of SIMNRA. The simulated areal 
density of the Ni film is 2.21 x 10!8 atoms/cm?, correspond- 
ing to a thickness of 242 nm. From the measured spectrum, 
the Ni peak area is A = 510000cts, and the areal densi- 
ty Nt = 2.20 x 10!8 atoms/cm’, corresponding to a thick- 
ness of 240nm, using the equation of Nt = A/oQQ, where 
the scattering cross section of Ni being o = 1.25 bar from 
do-/dQ = [Z| Ze? /(4E)|2(4/ sin? 8), and the number of inci- 
dent ions being Q = 4.7 x 10'9(7.5 uC). 
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Fig. 4. (Color online) RBS spectrum of a Ni film on SiO,/Si sub- 
strate, and the simulated spectrum using SIMNRA. 


C. Composition analysis of multi-element and multilayer 
samples 


A MoC thin film was prepared by multi-arc ion plating 
on Si substrate with Mo film as a transition layer. The Mo 
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Fig. 5. (Color online) The measured and simulated RBS-spectra 
(a) for 1.75 MeV protons backscattered from MoC film deposited 
on Si substrates with Mo as a transition layer, and RBS spectra of 
MoC/Mo/Si at 1.75-1.79 MeV(b). 


and C contents, and thickness of the MoC film, were deter- 
mined by RBS analysis using a 1.75 MeV proton beam and 
the backscattering angle of 9 = 170°. The RBS spectra are 
shown in Fig. 5(a) together with a fit using the SIMNRA code. 
The Mo signal is approximately a rectangular peak with the 
front edge rising sharply. It indicates that molybdenum in the 
MoC layer is quite uniformly distributed and the hump at the 
back edge is the signal of Mo from transition layer. The car- 
bon signal is a sharp peak rather than a hillock. Because non- 
Rutherford elastic backscattering occurred for carbon and the 
non-Rutherford backscattering cross section is 60 times larger 
than that of Rutherford backscattering at 1.74 MeV. [7] With 
1.75 MeV protons, non-Rutherford elastic backscattering of 
carbon occurred at depth near the surface. A raised Si edge 
can also be observed due to the enhanced elastic scattering 
cross-section for Si atoms at 1.75 MeV of protons. According 
to the simulation result, the content of Mo is 0.71 and the 
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Fig. 6. (Color online) Light elements analysis of TiBN film by non- 
Rutherford resonance scattering. 


content of C is 0.29 in the compound film layer with an area 
density of 1.52 x 10! atom/cm?. The Mo transition layer is 
296.4nm in thickness. Fig. 5(b) shows RBS spectra of the 
sample measured under energies of 1.75—1.79 MeV. It can be 
seen that the spectrum broadened and the C peaks moved to 
the low energy region with increasing incident beam energy, 
due to that occurrence of non-Rutherford elastic backscatter- 
ing of carbon deepened with the proton energy, and the energy 
resolution worsened with increasing proton energy. 


D. The detection and analysis of light elements in compound 
film 


It is quite difficult to detect films containing light elements 
on a substrate of higher atomic mass by backscattering tech- 
niques, in which the yield of substrate is too high to distin- 
guish signals of the light elements. In this case, one can use 
resonant scattering, which is of greater scattering cross sec- 
tion than Rutherford scattering. Appropriate higher energy 
should be applied so that resonance nuclear reaction occurs 
for certain elements. 

A sample of TiBN film deposited on Si substrate by multi- 
arc ion plating method was analyzed by 2 MeV proton with a 
backscattering angle of 0 = 170°. Fig. 6 is the energy spec- 
tra of the sample, with the fitting line and spectrum of the 
TiBN film after subtraction of Si substrate. The strong non- 
Rutherford resonances in 'H—'“N and 'H—"'B cross sections 
lead to sufficient enhancement of the N and B signals and 
areal density of the film can be calculated with the peak ar- 
eas. The net integrated counts of the Ti, N and B peaks in 
the blue spectrum are Ar; = 186921 cts, Ay = 118432 cts 
and Ag = 47248 cts, respectively. A large uncertainty may 
be caused by subtraction of the substantial signal due to the 
Si substrate. Accuracy of the N and B yields is determined by 
the error of intercept and slope of the straight regression line 
y = yo + kx, where y is an approximation of the yield Y in 
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channel x, yo is the intercept and k is the slope, best fitted into 
the experimental data. With 2 MeV protons, the cross section 
of Ti is Rutherford and the cross sections of N and B are non- 
Rutherford. The 'H—!4N cross-section enhancement factor is 
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(o/oR)n = 5 and the 'H—!'B factor is (7/op)y = 7 [7]. The 
cross section of Rutherford backscattering of Ti, B and N can 
be calculated as 


oR Ti(Eo, 170°) = (0.6359/27) x 107% = 1.590 x 107% cm?/sr, (1) 
oRN(Eo, 170°) = (0.0638/27) x 107% = 1.595 x 107% cm*/sr, (2) 
oRB(Ep, 170°) = (0.0323/27) x 10° = 8.078 x 10-7” cm?/sr. (3) 


The number of backscattered particles detected using the 
surface energy approximation for normal incidence of beam 
can be expressed as Ay = 0(E,)QQNt, where o(Eo) is scat- 
tering cross-section at incident proton energy Eo, Q is solid 


angle, Q is incident proton flux, N is number density of target 
atoms, and ¢ is film thickness. Then, for Ti: N and Ti: B ratios 
from the backscattered spectrum, we have 


Nri/Nn = (Ati/An) [oR N(Eo, 170°)/oR ti(Eo, 170°) (7 /or)N, (4) 


Nri/Ng = (Ati/As) [or B(Eo, 170°)/oRti(Eo, 170°) (7 /or)B- (5) 


The average stoichiometric ratio for this film was calculat- 
ed as Ti:N:B=1.4:1.75:1. 


IV. CONCLUSION 


Measurements were carried out by ion beam analy- 
sis establishment on the 1.7MV tandem accelerator at 
Wuhan University. The RBS spectrum calibration using a 
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